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Abstract. We review our theoretical work on the dynamics of a localized electron 
spin interacting with an environment of nuclear spins. Our perturbative calculation 
is valid for arbitrary polarization p of the nuclear spin system and arbitrary nuclear 
spin J in a sufficiently large magnetic field. In general, the electron spin shows rich 
dynamics, described by a sum of contributions with exponential decay, nonexpo- 
nential decay, and undamped oscillations. We have found an abrupt crossover in 
the long-time spin dynamics at a critical shape and dimensionality of the electron 
envelope wave function. We conclude with a discussion of our proposed scheme to 
measure the relevant dynamics using a standard spin-echo technique. 
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1. Introduction 

Experiments with trapped ions[l] and nuclear magnetic resonance (NMR) [2] 
have proven that the basic elements of a quantum computer can be real- 
ized in practice. To progress beyond proof-of-principle experiments, the 
next generation of quantum information processors must overcome sig- 
nificant obstacles regarding scalability and decoherence. The scalability 
issue is largely solved by proposals for a solid-state implementation 
for quantum computing, where established fabrication techniques can 
be used to multiply the qubits and interface them with existing elec- 
tronic devices. Due to their relative isolation from the surrounding 
environment, single electron spins in semiconductor quantum dots are 
expected to be exceptionally robust against decoherence. [3] Indeed, 
longitudinal relaxation times in these systems have been measured to 
be T\ = 0.85 ms in a magnetic field of 8T[4] and in GaAs quantum 
wells, the transverse dephasing time T| for an ensemble of electron 
spins (which typically provides a lower bound for the intrinsic decoher- 
ence time T2 of an isolated spin) has been measured to be in excess of 
100 ns. [5] 

Possible sources of decoherence for a single electron spin confined 
to a quantum dot are spin-orbit coupling and the contact hyperfine 
interaction with the surrounding nuclear spins. [6] The relaxation rate 
due to spin-orbit coupling 1/Ti is suppressed for localized electrons 
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at low temperatures [7] and recent work has shown that T2, due to 
spin-orbit coupling, can be as long as T\ under realistic conditions. [8] 
However, since spin-carrying isotopes are common in the semiconductor 
industry, the contact hyperfine interaction (in contrast to the spin-orbit 
interaction) is likely an unavoidable source of decoherence, which does 
not vanish with decreasing temperature or carefully chosen quantum 
dot geometry. [9] 

In the last few years, a great deal of effort has been focused on a 
theoretical description of interesting effects arising from the contact 
hyperfine interaction for a localized electron. The predicted effects in- 
clude a dramatic variation of T\ with gate voltage in a quantum dot 
near the Coulomb blockade peaks or valleys, [10] all-optical polarization 
of the nuclear spins, [11] use of the nuclear spin system as a quantum 
memory,[12, 13] and several studies of spin dynamics. [14, 15, 16, 17, 
18, 19, 20, 21] Here, our system of interest is an electron confined to 
a single GaAs quantum dot, but this work applies quite generally to 
other systems, such as electrons trapped at shallow donor impurities in 
Si:P.[9] 

An exact solution for the electron spin dynamics has been found 
in the special case of a fully polarized initial state of the nuclear spin 
system. [14, 22] This solution shows that the electron spin only decays 
by a fraction oc jj of its initial value, where N is the number of nuclear 
spins within the extent of the electron wave function. The decaying 
fraction was shown to have a nonexponential tail for long times, which 
suggests non-Markovian (history dependent) behavior. For an initial 
nuclear spin configuration that is not fully polarized, no exact solution 
is available and standard time-dependent perturbation theory fails. [14] 
Subsequent exact diagonalization studies on small spin systems [23] have 
shown that the electron spin dynamics are highly dependent on the type 
of initial nuclear spin configuration. The unusual (nonexponential) form 
of decay, and the fraction of the electron spin that undergoes decay may 
be of interest in quantum error correction (QEC) since QEC schemes 
typically assume exponential decay to zero. 

In the following section we describe a systematic perturbative theory 
of electron spin dynamics under the action of the Fermi contact hyper- 
fine interaction. Further details of this work can be found in reference 
[24]. 



2. Model and perturbative expansion 

We investigate electron spin dynamics at times shorter than the nu- 
clear dipole-dipole correlation time Tdd (?dd ~ 10~ 4 s in GaAs is given 



hfi2004coishloss.tex; 2/02/2008; 11:43; p. 2 



3 



directly by the inverse width of the nuclear magnetic resonance (NMR) 
line[25]). At these time scales, the relevant Hamiltonian for a descrip- 
tion of the electron and nuclear spin dynamics is that for the Fermi 
contact hyperfine interaction: 



where S is the electron spin operator, b gives the electron Zeeman 
splitting and h = J2k Ak\k is the quantum nuclear field. S± = S x ± iS y 
and h± = h x ± ih y are defined in the usual way. The hyperfine coupling 
constants are A^ = Avo\i^{rk)\ 2 where vq is the volume of a crystal unit 
cell containing one nuclear spin, tjj(r) is the electron envelope wave 
function, and A is the strength of the hyperfine coupling. While the 
total number of nuclear spins in the system, iVt t, may be very large, 
there are fewer spins (N < A^ot) within the extent of the electron 
wave function. In GaAs, all naturally occurring isotopes carry spin 
/ = |. In bulk GaAs, A has been estimated[25] to be A = 90 \ieV . This 
estimate is based on an average over the hyperfine coupling constants 
for the three nuclear isotopes 69 Ga, Ga, and 75 As, weighted by their 
relative abundances. Natural silicon contains 4.7% 29 Si, which carries 
I = \i and 95% 28 Si, with 1 = 0. An electron bound to a phosphorus 
donor impurity in natural Si:P interacts with N ~ 10 2 surrounding 29 Si 
nuclear spins, in which case the hyperfine coupling constant is on the 
order of A » 0.1 fieV. [9] 

For large magnetic fields 6, the flip-flop processes due to Vg are 
suppressed by the electron Zeeman splitting, and it is reasonable to take 
H ss Hq (zeroth order in Vg). The zeroth-order problem is algebraically 
simple, and leads to some interesting conclusions regarding the initial 
conditions. In this limit the longitudinal spin is time-independent, since 
[S z , Ho] = 0, but the transverse spin can undergo nontrivial evolution. 
Assuming uniform hyperfine coupling constants (Ak = A/N) and nu- 
clear spin I = 1/2 we evaluate the transverse electron spin dynamics for 
a nuclear spin bath of polarization p along the z-axis and two different 
nuclear spin initial states. When the Ak are uniform, the transverse 
spin exhibits periodic dynamics. However, at times much less than the 
period, given by the inverse level spacing (t <C N/A, setting h = 1), we 



H = H + V s , 

H = (b + h z )S z ; V s = ^(S + h_ + S-h + ), 



(1) 
(2) 



find 



in \unprcp. 
\ b +)t 




(3) 



(S + ) p t 



<S + > e" 



,iui n t 



co n = b + pA/2. 



(4) 
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In (3), the nuclear spin system is "unprepared". This state corresponds 
to either a translationally invariant direct-product state of the form 
life (V7l lt)fc + V 1 - /? U)fc)> wh ere p = 2/ T + 1, or to a statistical 
mixture of product states of the form IfjtT • • •), with average polar- 
ization p. In contrast, (4) corresponds to the case when the nuclear 
system has been "prepared" in an eigenstate \n) of the operator h z : 
h z \n) = [h z \ nn \n),(\h z \ nn = pA/2). It is important to note that the 
decay of the "unprepared" state is not irreversible. This decay can 
be recovered by performing a spin-echo measurement on the electron 
spin, since a 7r-rotation of the electron spin S z — > —S z reverses the 
sign of Hq: Hq — ► — Hq, and results in time-reversed dynamics. When 
the hyperfine coupling constants are allowed to vary in space, higher- 
order corrections in will, however, result in irreversible decay. In the 
following, we generalize to nonuniform and arbitrary nuclear spin / 
to evaluate this decay for an initial nuclear state that is an eigenstate 
of the /i^-operator. 

An exact generalized master equation (GME) can be derived for 
the electron spin operators beginning from the von Neumann equation 
for the full density operator (p = — i [H, p]).[26] The resulting GME 
is expanded in terms of VJf through direct application of the Dyson 
identity. We find, quite remarkably, that the equations of motion for 
the longitudinal ((S z ) t ) and transverse ((S+) t = (S x ) t + i {S y ) t ) spin 
are decoupled to all orders in Vg and take the form 

(S z ) t = N z (t) - i f dt'X zz (t - (S z ) t , , (5) 
Jo 

(S+) t = iu n (S+) t - i f dt>Z ++ (t - (S + ) t , . (6) 

Jo 

These integro-differential equations can be converted to a pair of alge- 
braic equations by Laplace transformation f(s) = J °° dse~ st f(t), Re[s] > 0. 
The algebraic equations yield 

S z (s) = i^o+^(£). s+{s) = (S + ) (?) 

S + l±j zz (S) S - lUJ n + I2u ++ {S) 

The numerator term N z (s) and self-energy S Z2 (s) are related to the 
self-energy matrix elements for spin-up and spin-down by N z (s) = 
-£(E TT (a) + E n (s)) and S 2Z (s) = S TT (s) - E n (a). All self-energy 
matrix elements are written in terms of an expansion in powers of 

V G : X a p(s) = E^(a) + sg( S ) + ---,a,P= (+,!,!)• For a sufficiently 
large Zeeman splitting b we find that all higher-order self-energy matrix 
elements are suppressed by a smallness parameter A: 

yMk+l)) () A fe. A — A m 

Kf3 (S)OCA , A- 2{b+pIA y (8) 
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In Born approximation (second order in the flip-flop terms Vg), and for 
an initial nuclear spin system with uniform polarization, we find 

(2) 

£jj(s) w £j-f(s) = — iNc+ [I+{s — ioo n ) + I-(s + iu n )\ , (9) 

E n (s) « £{f (a) = iiVc_ [7_(a - ««;„) + J+(s + iu n )] , (10) 

£ ++ (s) « E<g_( s ) = -iN [c_/ + (a) + c + /_(a)] , (11) 

c ± = /(/ + 1) - J- ]T <JJ(/jE ± 1)) , (12) 
A/tot fc 

J±(a) = — > K- « — / dx— ! r 1 -, f = 1.(13) 

In evaluating (13), we have assumed the electron is in its orbital ground 
state, described by an isotropic envelope wave function of the form 

r , / \ mi 

ip(rk) = V'(O) ex P ~2\l k ) • ^ ne index k gives the number of nuclear 



spins within radius r& of the origin and N is the number of nuclear spins 
within radius Iq, so in d dimensions (jjf) = 77- From the definition 

r / \ — " 

Ak oc \ifj(rk)\ this gives the coupling constants Ak = 2exp — f-^J d 

in units where A)/2 = 1. For times t < iV 3 / 2 /^, it is strictly valid to 
take iVt t — > 00 and change the sum to an integral in (13), [24] which 
gives the final result for I±(s), above. 

The time-dependent spin expectation values can now be recovered 
from the Laplace transform expressions by evaluating the Bromwich 
inversion integral: 

<*>< " h L " ae " SAsy - <s+) < = s? L <ise * s+(s) ' (14) 

To evaluate these integrals we close the contour in the negative real 
half-plane, avoiding all branch cuts of the functions S z (s), S+(a), and 
apply the residue theorem. The relevant contour is shown in Figure 1 for 
S+(s) or S z (s — iuj n ) within Born approximation when d = m = 2 (this 
applies to a two-dimensional quantum dot with parabolic confinement). 
The spin expectation values that result from this procedure are the 
sum of contributions with undamped oscillations (from poles on the 
imaginary axis), exponential decay (from poles with finite negative real 
part) and nonexponential long-time tails (from branch cut integrals). 
Since the contributions from poles on the imaginary axis do not decay, 
the spin expectation values (in a suitable co-rotating frame) are given 
generically in terms of a constant piece, and a remainder with nontrivial 
dynamics: 

(S x ) t = (S x )o + Rx(t), X = +,z. (15) 
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Figure 1. Contour integral (left) that determines all contributions (right) to the 
spin expectation values, including undamped oscillations (right, top), exponential 
decay (right, middle) and nonexponential decay (right, bottom; solid line: numerical 
integration, dashed line: analytic asymptotic expression). All three contributions are 
added to obtain p](t) — 1/2 + {S z ) t (right, inset). For these calculations we have 
chosen d = m = 2, / = 1/2, the initial condition pi(0) = 0, and values of b and p 
that correspond to A = 10/11. 



We show the different contributions to R z (t) in Figure 1 for d = m = 2 
in the weakly perturbative regime (where A < 1). 

In the strongly perturbative regime (when A <C 1), the time depen- 
dence of Rx(t) is given exclusively in terms of the functions I±(t) = 
2^1 J<7 dse st I±{s) that appear in (13), above. There is an abrupt crossover 
in the long-time behavior of these functions at a critical value of d/m. 
For d/m < 2, the major contributions to I±(t 3> N/A) come from 
the upper limit of the integral in (13), corresponding to nuclear spins 
near the center of the quantum dot. This leads to a modulation of the 
spin dynamics at a frequency A/2N. When d/m > 2, the major contri- 
butions come from the lower limit, corresponding to nuclear spins far 
from the dot center. In this case, the long-time dynamics are smoothly 
varying, with no modulation: 

f (i\ d/m piij&t v _ J. , < i 
R x (t) ~ I ± (t » N/A) oc W f ™ • ' - m 1 < 1 . (16) 

The difference in these two cases should be visible in a spin-echo 
experiment that uses a Carr-Purcell-Meiboom-Gill (CPMG) spin-echo 
sequence: | — (r — ir x — r — ECHO — r — ir_ x — r — ECHO) rep eat- 
We consider the strongly perturbative limit (A <C 1), and to resolve 
the relevant dynamics, the time between 7r-pulses must satisfy r <C 
v^dd, where 5 = A 2 /N, and r^d is the nuclear spin dipolar correlation 
time. [24] Under these conditions the CPMG echo envelope decay rate 
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Figure 2. Longitudinal-spin decay rate of the CPMG echo envelope 
(1/T 1 M oc R z (2t) /2t) as a function of the free evolution time 2r between 
7r-pulses. We plot the results for an electron bound to a phosphorus donor, where 
N — 10 2 (top) and a two-dimensional GaAs quantum dot with N = 10 5 (bottom). 

as a function of r is determined exclusively by the remainder term 
according to l/T^ 1 oc R z (2t)/2t for the longitudinal component and 
1/T 2 M oc Re[-R + (2r)]/2T for the transverse components. We plot the 
CPMG decay rate as a function of 2r in Figure 2 for two systems of 
interest. For an electron trapped at a donor impurity in bulk silicon, 
d = 3 and the orbital wave function is exponential (m = 1). This 
corresponds to v = 2 in (16). In a two-dimensional GaAs quantum dot 
(d = 2) with parabolic confinement, the ground-state orbital electron 
wave function is a Gaussian (m = 2), which corresponds to u = in 
(16). 



3. Conclusions 

We have reviewed our theoretical description for the dynamics of a 
localized electron spin interacting with a nuclear spin environment. We 
have predicted a sharp crossover in the relevant dynamics at a critical 
value of the dimensionality and form of the electron envelope wave func- 
tion, and have described a standard method that could be used to reveal 
the relevant dynamics. We stress that the electron spin dynamics are in 
general very rich, described by contributions with exponential decay, 
nonexponential decay and undamped oscillations. Furthermore, this 
work may have profound implications for the future of spin-based solid- 
state quantum information processing and quantum error correction, 
where previous studies have assumed exponential decay to zero. 
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